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£H ', Abstract. MAGIC currently is the largest imaging atmospheric Cerenkov 

telescope world-wide. Since 2004, 7-ray emission from several active galactic 
nuclei in the GeV/TeV energy range has been detected, some of which were 
newly discovered as very-high energy 7-ray sources. The 7-rays are assumed 
■ to originate from particle acceleration processes in the AGN jets. We give an 

overview of the AGN observed and detected by MAGIC, discuss spectral and 
temporal properties of these and show physics implications of some selected 
i-C ■ observations. 
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1. Introduction 

The study of very high energy (VHE, E > 100 GeV) 7-ray emission from active 
galactic nuclei (AGN) is one of the major goals of ground-based 7-ray astronomy. 
ON \ The sensitivity of the current imaging air Cerenkov telescopes (IACT) enables 

phenomenological studies of the physics inside the relativistic jets in blazars, and 
in particular advances in understanding both the origin of the VHE 7-rays as 
well as the relations between photons of different energies (from radio to VHE). 

Except for t he radio galaxy M87, all 17 currently known VHE gamma-ray 
emitting AGNs (jWagnerl l2007al lbl) are blazars Q which are characterized by a 
close orientation of the jets with the line of sight. Their spectra are dominated 
by non-thermal emission that consists of two distinct broad components. While 
• '— j . the low energy bump, located at optical to X-ray energies, is unisonously ex- 

plained by synchrotron emission of electrons, the origin of the high-energy part 
of the spectrum is still debated. Leptonic models ascribe it to inverse Comp- 
ton processes that either up-scatter synchrotron photons, or to external photons 
that originate from the accretion disk, from nearby massive stars, or are re- 
flected into the jet by surrounding material. In hadronic models, interactions of 
a highly relativistic jet outflow with ambient matter, proton-induced cascades, 
or synchrotron radiation off protons, are responsible for the high energy photons. 

Another defining property of blazars is the high variability of their emission 
ranging from radio to 7-rays. For VHE 7-ray blazars, correlations between X- 
ray and 7-ray emission have been found o n time scales ranging from several 10 
minutes to days (e.g., iFossati et al. 20041 ) . although the relation has proven to 
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be rather complicated. 

1 See http://www.mppmu.mpg.de/~rwagner/sources/ for an up-to-date source list. 
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Here we present selected results for MAGIC blazar observations of Mkn 421, 
Mkn 501, 1ES 2344+514, PG 1553+113, and BL Lacertse. 



2. The MAGIC Telescope 



MAGIC (jBaixeras et al.l l2004i ; ICortina et"all 120051 ) is located on the Canary 
Island of La Palma (2200 m a.s.l.) and is currently the largest IACT. Its energy 
range spans from 60 GeV (trigger threshold at small zenith angles) up to tens 
of TeV. MAGIC has a sensitivity of ~ 2.5% of the Crab Nebula flux in 50 
observation hours. Its energy resolution is about 30% above 100 GeV and about 
25% fr om 200 GeV onwards . The MAGIC standard analysis chain is described 
e.g. bv lWagner et al.l ( 20051 ). Obser vations during mode rate moonshine enable 
a substantially extended duty cycle ( Albert et al.l 2007d ). which is particularly 
important for blazar observations. Parallel optical observations are performed 
with the KVA 35 cm telescope. 

The MAGIC AGN observation program encompasses over 5 00 hours /year. 



Havashida et al 



Mazin et al.l 



Besid es participating in various multiwavelength campaigns (e.g. 
2007 ) and many target of opp ortunity observation agreements (e.g., 
2007 : lAckermann et al. 20071 ). MAGIC is cooperating with other Cerenkov tele- 
scopes a s to extend the dy namical range and/or the time coverage of obser- 



vations ( Mazin et al. 2005) and monitors some 

^ 120071 1. 



a regular basis (jGoebel et al. 



well-known bright blazars on 
Apart from having observed previously 
known blazars, MAGIC has also newly disc overed a co uple of blazars in VH E 7- 
rays, as 1E S 1218+30.4 (I Albert et al.l l2006aTl. Mkn 180 (lAlbert et al.ll20q6bl'l. P G 
1553+113 (|Albert et al.ll2007ah . and recently BL Lacertse (jAlbert et al.ll2007fh . 



3. Markarian 421 

Mkn 421 is the closest (z = 0.030) known TeV blazar and, along with Mkn 501 
also the best studied one. It has shown variations larger than one order of mag- 



nitud e and occasional flux doubling times as short as 15 minutes (IGaidos et al 
119961 ). Variations in the hardness of the TeV 7-ray spe c trum during flares wer e 
reported by several groups (e.g., iKrennrich et al.ll2002l : lAh aroman et all 120051 ) 



Simultaneous observations in the X-ray and VHE bands s howed a significant 
correlation of the respective fluxes ( Krawczvnski et al.ll200l1 ). 



Mkn 421 was observed with MAGIC during 19 nights for a total of 25.6 
hours, with observation times per night ranging from 30 minutes up to 4 hours. 
It was found to be in a medium flux state ranging from 0.5 to 2 Crab units 
above 200 GeV. Significant variations of up to an overall factor of four and 
up to a factor two between successive nights were found. A clear correlation 
between the 2 — 10 keV X-ray count rate, obtained from the public RXTE All- 
Sky Monitor data0 and the VHE 7-ray flux can be described well by both a 
linear fit (correlation coefficient r = 0.641q22) and a parabolic fit (both forced 
to go through the origin). 



2 These data are available at http://xte.mit.edu/. 
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Figure 1. The measured 
(grey squares) and intrin- 
sic (black points) spectrum 
of Mkn 421. Solid line: 
power-law + cutoff fit to 
the intrinsic spectrum. 



Figure 2. Correlation be- 
tween spectral slope and 
differential flux for various 
VHE observations at f TeV for 
Mkn 42 f (for ref erences, see 
lAlbert et al.ll27)07ch . 



The energy spectrum of Mkn 421, for the first time measured down to 
E = 100 GeV, is shown in Fig.[T]along with the reconstructed intrinsic spectrum. 
The latter has been obtained by correcting for 77 absorption effects caused by 
extragalactic background ligh t (EB L; Hauser & Dwek 2001), using the recent 
EBL model of Primack et al.l ( 20051 ). The intrinsic spectrum is clearly curved, 
implying that the curvature in the measured spectrum has an intrinsic origin 
rather than being caused by EBL absorption. Both a power-law+exponential 
cutoff fit and a curved power-law fit indicate a flattening of the spectrum towards 
100 GeV. The spectrum is hard ening with an i ncreas ing flux level (Fig. [2]). The 



analysis is described in detail in lAlbert et al.l ( 2007d ) . 



4. Markarian 501 

Mkn 501 (z = 0.034) is known to be a strong and variable VHE 7-ray emitter. 
During a flare in 1997 Mkn 501 showed strong variability on timescales of 0.5 
days, and th e integral flux reached 10 times the flux of the Crab nebula above 
1 TeV (e.g., Aharonian et al. 19991 ). While a correlation between spectral hard- 



ness and 7-ray flux was found, the position of the inverse Compton peak could 
not be detected. 

MAGIC observed Mkn 501 for 24 nights during six weeks in summer 2005. 
About 18 hours out of the total 30 observation hours were performed in the 
presence of (moderate) moonshine. During most of the observations the source 
was in a rather low flux state, corresponding to 30% of the Crab nebula flux at 
E > 200 GeV. In two nights (one with moon present), however, Mkn 501 was 
reaching up to four times the Crab nebula flux. A correlation between spectral 
hardness and VHE 7-ray flux was found also for Mkn 501 (Fig. [3|). For the flare 
nights, a clear peak in the spectrum is seen (e.g., Fig. 0]). In leptonic models it 
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Figure 3. Correlation of spectral hard- 
ness and VHE 7-ray flux above 200 GeV. 
Each data point represents the average 
over one night, except for the four high- 
est VHE flux points (two points per 
night). 



Figure 4. Measured 
(lower curve) and recon- 
structed intrinsic (upper 
curve) spectrum of Mkn 
501 during the night of 
2005 June 30. 



can be identified with the "Inverse Compton peak" , which now for the first time 
is seen unambiguously in VHE range for blazars. 

A closer look at the two flare nights reveals rapid flux changes with dou- 
bling times as short as 3 minutes or less (Fig. [5]), while constant background 
rates assure that no instrumental effects are responsible for these. Fig. [6] shows 
corresponding light curves in different, distinct energy bands between 150 GeV 
and 10 TeV. Interestingly the flares in the two nights behave differently: While 
the 2005 June 30 flare is only visible in 250 GeV to 1.2 TeV, the 2005 July 09 
flare is apparent in all energy bands, and a marginal time delay towards higher 
energies is visible. For the first time, short (~ 20 min) flares with a resolved 
time structure can now be subjected to detailed studies of particle acceleration 
and coo ling timescales. The analysis and more results are shown in I Albert et al. 
(j2007dh . 



5. 1ES 2344+514 

1ES 2 344+514 (z = 0.044) was detected in VHE 7-rays in 1995/6 on the 5.8 a 
level (ICatanese et all 1 19981 ). Over 1/3 of the excess was recorded during one 



night in December 1995 with F 7 (> 350GeV) = 64% of the Crab nebula flux. 
The 7 -ray emission was confirmed on a 4.4cr significance level bv lAharonian et al 
(120041 ) in observations in 1997-2002. 



In 2005, MAGIC saw a clear lla excess from the direction of 1ES 2344+514. 
The measured flux corresponds to 10% of the Crab nebula flux {E > 200 GeV). 
For the first time a diurnal measurement of the emission of this blazar was possi- 
ble, yielding no strong evidence for significant variability on times scales of days 
and thus giving first insights to the day-scale properties of a low blazar emission 
state. The differential energy spectrum can be fitted by a simple power law with 
a photon index of a = —2.95 + 0.12. The derived spectrum is s ofter than the one 



repor ted by the Whipple collaboration duri ng the flare in 1995 (ISchroedter et al 
2005). Details on the analysis are given in Albert et al. ( 2007bl V 
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Figure 5. Light curves (2 min bin width) for the nights of 2005 June 30 
(top left) and 2005 July 9 (top right). The lower plots show the corresponding 
background rates, which were found to be constant. 
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Figure 6. Evolution of the flares of 2005 June 30 (left column) and 2005 
July 9 (right column) in distinct energy bins. 
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Figure 7. E > 180 GeV 
and i?-band light curves 
for PG 1553+113. 



I™ 4 

LU 

£io- 11 



PG 1553+113, MAGIC 2C05+2C06 
PG 1553+113, H.E.S.S. 2005 
Crab Nebula, MAGIC 2004+2005 




dF/dE = f„- (E/200 GeV) 



P f = (0.18 ±0.03) ■ 10 cm- 2 s- 1 TeV 
T = (4.21 ± 0.25) 



500 E [GeV] 



Figure 8. Differential energy 
spectrum for PG 1553+113. 



6. PG 1553+113 

Based on its SED properties, PG 1553+113 was one of the mo st promising 
candidates for VHE 7-ray emission ( Costamente &: Ghisellini 20021 ). The blazar 
was observed in 2005/6, yielding a 7-ray signal at a significance level of 8.8 a. 
There is no evidence for short term variability on a time scale of days, but 
a significant change by a factor of three in the flux level from 2005 to 2006 
was found. The VHE 7-ray flux was not found to correlate with a substantial 
brightening of the source in the optical i?-band around MJD 53820 (Fig. [7J). 
The differential energy spectrum for PG 1553+113 is well described by a pure 
power law with a photon index of a = —4.2 ± 0.4, the steepest photon index 
observed so far for blazars (Fig. [SJ. The signal detected by MAGIC confirms a 
tentative signal seen by H.E.S.S. at a high er energy threshold wi th data taken 
at about the same time as MAGIC in 2005 (jAharonian et al.ll2006ft . The source, 
therefore, can now be considered as a firmly e s tablish ed VHE 7-ray emitter. 
Details of the analysis are given in lAlbert et all l|2007d ). 

The distance of PG 1553+113 is currently unknown, despi te several at- 
temp ts of determining its redshift (Direct measurements: z > 0.09, ISbarufatti et al 
2006j). The objects shows no emission lines, which probably implies a close align- 
ment of the jet with the line of sight, or a jet outshining the galaxy core. Using 
the observed spectrum of PG 1553+113 a nd assuming a minim um possible den- 
sity of the evolving EBL (lower limit in Kneiske et al. 20041 ). an upper limit 
of z < 0.74 can be inferred wh en requiring the intrin s ic photon index T < 1.5. 
Other works employ SED fits (lAharonian et al.ll2006l; I Albert et al.ll2007ah . the 
structure o f the intrinsic s pectrum ( Mazin &: Goebelll2007l ). or blazar ensemble 
arguments (|Wagnerl l2007al ) to determine redshift limits. 



7. BL Lacertae 

BLLacertae (z = 0.069) is the historical prototype of the BL Lac class of AGNs. 
Unlike for the up to recently known "VHE blazars" with their synchrotron peaks 
in the X-ray domain, the peak of its synchrotron emission is located in the sub- 
millimeter to optical band (^pggjj = 2.2 x 10 14 Hz). 
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22.2 hours of MAGIC observations of BL Lacertas from August to December 
2005 resulted in a 5.1a VHE 7-ray signal. Above 200 GeV, an integral flux of 
(0.6 ±0.2) x 10~ n cm -2 s _1 was measured, corresponding to approximately 3% 
of the Crab nebula flux. No significant evidence of flux variability was found. 
Optical and radio data shows enhanced flux levels during the observations; While 
the VHE 7-ray flux shows similar variations, no correlation can be claimed at 
this point. 

The differential spectrum between 150 and 900 GeV was found to be rather 
steep, with a photon index of a = —3.6 ± 0.5. In contrast to a 1997 flare of 
BLLacert ae, which necessitat ed external photon fields to explain the VHE 7 
emission (Ra vasio et al.l 20021 ) , the MAGIC observati ons can be ac c ommo dated 



in a simple SSC model. The detection is described in lAlbert et al.l (20073). 
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